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We study the Sivers single spin asymmetry of the Drell-Yan lepton pair production in the small- a; 
regime. We find that in the corresponding kinematic region the spin asymmetry calculated in the 
small- j: approach is consistent with either the usual transverse-momentum-dependent factorization 
formalism or the coUinear factorization formalism, respectively. We estimate the Sivers asymmetry 
for both polarized p-f p and p-f A collisions and argue that the Drell-Yan production is an interesting 
. . . and unique probe for both the transverse spin physics and the small-x saturation effect. 
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O ' I. INTRODUCTION 

D 

. Physics at forward rapidities at Relativistic Heavy-Ion Collider (RHIC) has attracted great attention from both 
0^ ' experimental and theoretical sides in recent years. It provides great opportunities to study novel hadronic physics 
phenomena, among which the single transverse spin asymmetry (SSA) in polarized p-l-p collisions [l[, and the small- a; 
gluon saturation in p-l-A (d-t-A) collisions are particularly interesting. On one hand, it has been realized that the 
<~i ' SSAs observed in high energy collisions are directly connected to the transverse motion of partons inside the polarized 
hadron. Thus it could help map out the three-dimensional image of the hadron in the transverse momentum space Q . 
^ ^ On the other hand, the single hadron suppression and the dihadron correlation in the forward d-|-Au collisions have 
^ , provided a way to study the small-x gluon saturation or multiple parton scattering effect 0-0. An effective theory 
' of the color-glass-condensate (CGC) has been applied to describe these phenomena. Seemingly completely different 
research topics, the SSAs and the small- a; gluon saturation can be closely related to each other and be studied in the 
I ' same process ^j. 

J> \ The SSAs have been observed in many experiments at different collision energies. It has been found that the SSAs 
^\ • become the largest in the forward direction of the polarized hadron [l| , where the incoming partons from the polarized 
hadron have relatively large momentum fraction x, thus mostly valence quarks. At the same time, the partons in the 
target have very small x in high energy collisions, and thus are dominated by the gluons and could be described by 
' the small- a; gluon saturation physics. In other words, the polarized p^-l-p and p^-|-A in the forward rapidity region at 
CNJ high energy collisions will be a unique opportunity to probe both the transverse spin and small-x saturation physics. 

The Drell-Yan production in high energy polarized p^-l-p and p^-f A collisions will be very important in such a 
, direction. The Drell-Yan production at small transverse momentum can be described by a well-established transverse 
momentum dependent (TMD) factorization formalism [lol - [l3j and thus the experimental data can be clearly analyzed 
and be cleanly interpreted. From the small- a: perspective, because of the lack of both final state interactions and 
fragmentation effects, Drell-Yan lepton pair production in p-l-A collisions is a very clean probe for the unintegrated 
gluon distribution, particularly the so-called dipole gluon distribution at small- a;. From the spin physics perspective, 
. Drell-Yan lepton pair production at small transverse momentum of the pair provides a way to study the so-called 
transverse momentum dependent distributions (TMDs) in the polarized hadron, encoding important information 
about the hadron structure beyond what has been learned from the usual collinear parton distribution functions 
(PDFs) . The particularly interesting one is the quark Sivers function [T^ , which represents a distribution of unpolar- 
ized quarks in a transversely polarized nucleon, through a correlation between the quark's transverse momentum and 
the nucleon polarization vector. Unlike the collinear PDFs which are universal, the TMDs or the Sivers functions are 
not exactly universal [T5| . instead, they have the time-reversal modified universality. It was shown from the parity and 
time-reversal invariance of QCD that the quark Sivers function in semi-inclusive deep inelastic scattering (SIDIS) and 
that in Drell-Yan process are exactly opposite to each other [l3| . The sign change of the Sivers function between 
the SIDIS and Drell-Yan production is a unique prediction of our current QCD TMD factorization formalism, and 
provides a critical test of the TMD factorization approach and our understanding of the SSAs flS]. 

The quark Sivers functions have been measured in SIDIS process at HERMES, COMPASS, and JLab experiments 
[l9l - l2l| . Future measurements of the SSAs in Drell-Yan production have been planned [S^-^ll, in the hope of verifying 
the sign change in the near future. In this paper, we study the Sivers asymmetry of the Drell-Yan production in high 
energy collisions, particularly in the small- a; regime. We first formulate both the spin-averaged and spin-dependent 
cross section in the small- a; formalism or CGC framework in Sec. II. We then compare our small-x formalism with 
those obtained in the TMD factorization in Sec. HI. We find that the two formalisms are consistent with each other 
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in the so-called geometric scaling region, where the invariant mass of the lepton pair M is much larger than the 
saturation scale Qs but saturation effects are still important. In Sec. IV, we further show that the small- a; formalism 
is also consistent with the usual coUinear factorization formalism in the so-called forward limit. In Sec. V, we present 
our numerical predictions for the Sivers asymmetry of the Drell-Yan lepton pair production in the forward rapidity 
region at RHIC energy for both p-|-p and p-|-A collisions. We summarize our paper in Sec. VI. 



II. SIVERS ASYMMETRY OF DRELL-YAN PRODUCTION IN SMALL-a; REGIME 



In this section, we will study the Drell-Yan lepton pair production in both transversely polarized p^+p and p^+A 
collisions, 



p\P,s^)+A{Pa)^[i*'^] £+£-{q)+X, 



(1) 



where P and s± is the momentum and the transverse spin vector of the polarized proton, A represents either a proton or 
a nucleus with a momentum Pa, and q is the momentum of the lepton pair with invariant mass = M"^ . To calculate 
the Drell-Yan cross section in small- a; regime, we first study the differential cross section in the high energy scattering 
of a quark off a color glass condensate, q{k)-\-A — ^ ^*{q)+X . The first step in the calculation is to obtain the light-front 
wave functions of the incoming quark splitting into a quark and a virtual photon, q(k, a) — > ^*{q, X) + q{k — q, fi), with 
a and /3 the helicity for the incoming and outgoing quarks, respectively. The splitting wave function in momentum 
space is given by 



<t>iiAk,q) 



1 



up{k - (7)7^e^(q, \)ua{k) 



^/8{k-q) + k+q+ {k-q)-+q- 



(2) 



where k^ = [fc+, fc^/2fc+, fc^] is the momentum of the incoming quark in the light-cone component, and q is the 
momentum of the virtual photon. For virtual photon with invariant mass Af, we could write the momentum q^^ as 



Ml 
2^ 



,q± 



with transverse mass M±_ = \/q^+ M^. We choose the polarization of the virtual photon as 



et^{q,X) 







e't(q,X) = — 



ql-M^ 



2q^ 



,q± 



(3) 



(4) 



where and are the transverse and longitudinal polarization vectors, and e^~^'^ = (=Fl, — i)/>/2. Our choice 
for the transverse polarization is consistent with those in Ref. [4^. With our choice for the longitudinal polarization 
vector, one could easily verify that the polarization vector satisfies the requirements: 



q ■ e{q, A) = 0, ^ e'^(g, A)e*''((Z, A) = -5'^"^ 



qf'q" 



(5) 



The cross section can usually be written as a compact form in the transverse coordinate space, in which the splitting 
wave function ^^p{k,q~^ ,r) is defined as the following: 

(6) 



V'^^(fc,<Z+,r) = J d'q^e^---<j,ip{k,q). 
With the polarization vectors defined above, we could easily calculate il}^g{k,q'^ ,r). The result is given by 



^J/(fc,(7+,r) = 27r^j^e'''^-'-ieMKi{eM\r\ 



^ [Sa-Sfj- + (1 - z)Sa + Sp+] , A = 1, 



^ [Sa+6p+ + {I - z)6a-6p-] , A = 2. 



(7) 



(8) 



where z = q^/k'^ and e\j = (1 — z)A'P, and KQ i{eM\r\) are the modified Bessel functions of the second kind. It 
is important to realize that we have kept the transverse momentum A:^-dependence for the incoming quark, which 
is crucial as the spin effect comes in as a quark TMD distribution inside the polarized proton. If we set k± = 0, 
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FIG. 1. An example diagram to illustrate the interaction between the quark and the target. Interactions before and after the 
splitting have to be taken into account for both the amplitude and the conjugate amplitude. 



these splitting wave functions reduce to those derived in Ref. |25| where a collinear incoming quark distribution was 
considered. 

Following Ref. [1, H^, when high-energy partons scatter off the CGC, the interaction is eikonal in that the 
projectile partons propagate through the target without changing their transverse position but picking up an eikonal 
phase. This phase can be expressed by the relevant Wilson lines. With the splitting wave function in the transverse 
coordinate space in hand, we could easily derive the differential cross section for a quark scattered off the CGC (see 
Fig. HI, q{k) + A '^*{q) + X, and the result is given by 



da{qA j*X) 
dq+(Pq±_ 



d^b d^r d^r' 



(27r)2 (27r)2 (27r)2 



(9) 



where is the quark fractional charge, v — zr + {1 — z)h and v' = zr' + (1 — z)h, and the two-point functions of form 
S^'^\x,y) take care of the multiple scattering between the 7*(7-pair and the nucleus target, and they are characterized 
by the Wilson lines with the following definition 



S^^l{x,v) = —{Tr{U{x)U\v)))^ 



(10) 



where the notation (• • • )xa is used for the CGC average of the color charges over the nuclear wave function and xa 
is the smallest fraction of longitudinal momentum probed, and is determined by the kinematics. U {x) is given by 



U(x) 



Vexpl^igsJ dx+T''A-{x+,x±)Y 



(11) 



with T'^ the generator of color SU{Nc) group in the fundamental representation, and V denoting an ordering in x'^ . 

Eq. ([9]) is a very compact form in the transverse coordinate space. In order to transparently incorporate spin- 
dependent effect and also to demonstrate the connection to the usual TMD factorization formalism for Drell-Yan 
production, we will transform Eq. (|9]) to the transverse momentum space, in which the cross section can be written 
in the following form 



da{qA j*X) 
dyd'^qjL 



= J d^ bd^p±_ F{xA,P±) [Ht iq±,kj_,pj_,z) + HL{q±,k^,pA_,z)], 



(12) 



where y is the rapidity of the virtual photon, and F{xa,P±) is the unintegrated gluon distribution (or dipole gluon 
distribution) defined as 



Fi^A,P^)= I S2^'^"-'^"]J^(Tr(f/(0)t/^(r^))), 



(13) 



Ht and represent the corresponding hard-part functions for transverse and longitudinal polarized virtual photon, 
respectively, and they are given by 

(14) 



HTiq±,k^,p^,z) - [l + (l-z)2] 
HL{q±,k^,p^, z)^ 2(1 -zfM^ 



q± - zk^ 



zk 



zpi 



{q±_ - zki_Y + . 



1 



_(g_L - zkj_Y + e 



if (^-L - zki_ - zp±_)^ + elj 
1 



(15) 
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To obtain the differential cross section for Drell-Yan production in polarized p^+p and p^+A collisions from above 
formalism of q + A — >■ 7* + X, one realizes that the unpolarized quark distribution inside the transversely polarized 
proton can be expanded as follows [28^ 



fq/pr{x,k±) = fq/p{x,k\) + ^-^;j^.ftT{x,k\), (16) 



p 



where Mp is the proton mass, and e"*^ is a two-dimensional anti-symmetric tensor with e^^ ~ 1. fq/p{x,k\) is the 

spin-averaged quark distribution function of flavor q, while f^j,^[x, k\) is the quark Sivers function. Thus by simply 
convoluting the spin- independent quark distribution function /^/^(a;, fc^) with Eq. (jl2p . we immediately obtain the 
spin-averaged cross section da = [(iCT(sj^) + (icr(— sj_)]/2 for virtual photon production, p^{P^ sx_) + A{Pa) ')*{q)-\-X, 

dyd^q^ ^2^^^" J. T ^^"^h/iM^^iJ J dbd p±F{xA,P±) 

X [HTiq±, k±,p±,z) + HLiq±, k±,p±,z)], (17) 

where runs over all light (anti)quark flavors, the quark momentum fraction x — Xp/z with Xp = M±/y/se^ , and 
the gluon momentum fraction in the target is given by xa = M±/ ^e~'^ . Eq. PT)) is consistent with the results 
derived before in [2^, if one sets the quark transverse momentum from the proton side k± — >■ 0. In our case, since 
we are interested in the spin-dependent effect which involves TMDs, we generalize the formalism to include the k±- 
dependence. Now convoluting the quark Sivers function with Eq. (|12p . one obtains the spin-dependent cross section 
dAa = [da{s±) — da{—s±)]/2 as 

J-p^^ = ^^'^ -d fc,^^.A/(x,fcJ j d bd p^FixA.P.) 

X [Ht iq±,kj_,p^,z) + HL{q_L,kj_,p^,z)]. (18) 

The differential cross section for Drell-Yan dilepton production (decayed from virtual photon) can be easily deduced 
using the following relation, 

dajpU £+i-X) _ aen. dajpU ^ 7*^) 

dM^dyd^qi_ iirNP dyd^q± ' ^ ' 

From Eqs. (IT71) . ([T5| and ([T5|) . we obtain the conventionally defined SSA (or the Sivers asymmetry) for Drell-Yan 
dilepton production as follows 

_ dA(7{p'^A^ £+£-X) / da{p^A^ £+£-X) 
^~ dAPdyd^q± / dAPdyd^q± ' ^ ' 

III. CONNECTION TO TMD FACTORIZATION FORMALISM 

In this section, we investigate the connection between the Drell-Yan differential cross section obtained above in the 
small- a: formalism and those from the TMD factorization approach. To find such a connection, we need to study the 
factorization property of the above differential cross section in the kinematic region where M 3> q_L. Furthermore, 
we assume that M is also much larger than the saturation scale Qs which sets the transverse momentum p± of the 
unintegrated gluon distribution F{xatP±)- In other words, we are going to study the region M ^ q± ^ k± ^ p± 
where the TMD factorization is supposed to hold for Drell-Yan production. 

In order to extract the leading power contribution from Eqs. (IT7|) and (US]), we notice that the integral in both 
equations are dominated by the end point contribution of z '-^ 1 where is in order of q^ [33,[3l|. Following Ref. [30| . 
we introduce a delta function: / dzS{z — 1/(1-1- A^/e|^)) = 1, where A^ = (1 — z){q± — k±)'^ + z{q± — k± — p±)'^ , and 
integrate out z first. This delta-function can be further expanded in the limit of q± <^ M, 

^l-z \ S{l^z) , S{l-z) ] 

z [(l-£)+ + (l-z)+J' ^'^^ 
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where a logarithmic term in the above expansion is power suppressed and has been neglected. Substituting this 
expansion to Eqs. ([T7)) and (|18p . we obtain the leading contribution to the differential cross section in the limit of 

q± < M, 



dM'^dy(Pqi_ 67r3M4 



j d^k^Xpfq/p{xp,k\) J dz J d'^bd^p±F{xA,P±)A{p±,q± - k±, z), (22) 



dAa{pU^i+i-X) V- 2 f ,2, e'^^ s^k'i 



dM^dyd'^qjL 6t:^M^ 



2 f 2 s"fc_L J- 2 



^ J dz J d^bd^pj_F{xA,Pj_)A{pj_,qj_- k±,z), (23) 
where A(pj^, ^j^, z) with = 9^ — A:^ is given by 

e±\i±-p±\ t^-Pi_ 



A{p^,t^,z) 



il^z)el + zi£^-p^y \£^-p^\ 



(24) 



Several comments are in order at this point. First, we found that the contribution from the second delta function 
i5(l — z) is power suppressed. To see this more clear, one can substitute e^^ = zA^/ (1 — z) for both hard-function Ht 
and Hl in Eqs. p7)) and ^TE\\ . there is an overall factor (1 — z)^. Second, we have also found that the contribution 
from hard-function (corresponding to longitudinal polarized virtual photon) is power suppressed compared to Ht 
(corresponding to transverse polarized virtual photon). So the final leading power contributions in Eqs. (|22l) and ([25)) 
are purely connected to the transverse polarized virtual photon. This is consistent with the expectation of the TMD 
factorization formalism, in which the virtual photon is generated from qq annihilation process and thus is transversely 
polarized. 

On the other hand, in the TMD factorization [lol - [l3| . we have the spin-averaged cross section as 



daip^A £+£-X) Ana, 



2 

cm 



dM^dyd'^q±_ ZNcM^ 



[ d^k^d^i^d^X^S"^ (fci +t^ + \^- q^) 



'XXpfq/p{Xp,kl^)xAfq/A{xAj±)H{M^,Xp,XA)S{X±), (25) 

and the spin-dependent cross section as 

dAaip'^A £+£-X) iiral^ 



H H / d^k^d^i±d^>^±S^ ik± +£± + X±- q±) 



dA'Pdyd'^q± iNcM'^ 

X '-^^^^P^xpf^^\xp, kl)xAfq/AixA, £l)H{M^, xp, xa)S{X^). (26) 

Here f^'^{xp,k'^), fq/p{xp,k'^) and fq/A{^Aj^'±) are the TMD quark Sivers function, quark and anti-quark distribu- 
tion, respectively. H{M'^,Xp,XA) and S{Xj_) are the hard and soft factors with the lowest order expressions: 

HiM^,Xp,XA) = l, SiX±)=S\X^). (27) 
To find the connection between the above TMD formalism and the small-x formalism, we can calculate the (anti)- 




FIG. 2. The Drell-Yan production in TMD factorization with the anti-quark distribution generated from the dipole gluon 
distribution. 
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quark distribution function fq/Ai^Aj^^) in terms of the dipole gluon distribution F{xa,P±) at small-x limit. This 
has been done in [30], and in the small-x Hmit, the result can be written as the following 

fg/A{xA,il) = ^J^I dH£p^F{xA,Pi_)A{p^,l^,z) (28) 

where z) is given by Eq. (|24p . At the same time, since now the TMD (anti)-quark distribution starts with 

nontrivial leading order expansion, we could set both the hard and soft factors as their corresponding leading order 
expressions in Eq. (|27p . Now feed the expression for the TMD quark distribution in Eq. ([^5]) back to the TMD cross 
sections in Eqs. (|25]) and (|26)) . taking the advantage of the delta function 5"^ {k± + i± + X± ~ q±) to integrate out 
(P£±, thus the factor A{p±,£j_, z) becomes A{p±, q± — k±, z) which also appeares in the small-x formalism in Eqs. ((22|) 
and (j23p . We immediately find that both the spin-averaged and spin-dependent cross sections agree completely with 
those in Eqs. (1221) and (|23p from the small-x formalism. For an intuitive understanding, see Fig. [2] where we show 
the Drell-Yan production in TMD factorization but with the anti-quark distribution generated from the dipole gluon 
distribution. Therefore, we have demonstrated that the small- a; calculation for the Drell-Yan differential cross section 
is consistent with the TMD factorization at this particular order. 



IV. CONNECTION TO COLLINEAR FACTORIZATION FORMALISM 



In this section, we discuss the connection between the differential cross section obtained in the small- a; formalism 
and those in the usual collinear factorization. For this purpose, we are interested in the kinematic region where 
M^, ^ ~ p\, k^. This is the region where the parton density in the nucleus target is still dilute and 
M^, q\ ^ Aqcd thus the usual collinear factorization is supposed to hold. 



A. Spin-averaged cross section 

We will start with the spin- averaged cross section. Since we are interested in the region where Af^, ^ fc^, thus 
we could drop the fc^ -dependence (compared to q±) in the hard-part functions Ht and Hl in Eq. p7)) . Once k± is 
dropped, the only k± dependence is then coming from the fq/p(x,k\), and we could integrate over cPk±^ to obtain 
the usual collinear parton distribution funciton, 

d?k^fq/p{x,kl)= f^/p{x). (29) 

On the other hand, for the nucleus target side, in the dilute parton region, we have [32j ] 

tfb(fp^p\F{xA,P±) = ^ " XAfg/AixA), (30) 

where fg/A{xA) is the usual collinear gluon distribution inside the nucleus. In order to obtain the needed p\ such 
that we could connect to the collinear gluon distribution from Eq. (|17p . we thus need to make the expansion with 
respect to p± in the hard-part functions Ht and Hl in the region where M^, q± ^ ^ p\, 



J (fb J (fp^F{xA,P±)p'[p'i^ g^p^Q^„ [i?T(g_L,fc_L = 0,P±,z) + HLiqA_,kj_ = 0,pj_,z)] 



(31) 



After the p± expansion, we set p± — > 0, thus the hard-part functions are now independent of p±. We could then use 
J (fb(fp±F{xA,p±)p'[^p'i^ J d'^bd'^p±F{xA,p±)p'i^g'^ = ^^^^^XAfg/A{xA)^g''^ . (32) 
Eventually we could write the spin-averaged cross section as 

dAPdyd^q^ =^^^K^Y^ j^--U/p{->Afg,A{xA)H{q^^zl (33) 

where the hard-part function H{q±, z) is given by 

I ^ d 

H{q±, z) = -g'^ [HT{q±, k± = 0,pj_, z) + HL{q±, k± = 0,p_L, z)]^^^^ (34) 
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With the Ht and given in Eqs. (1141) and (I15|) . we could immediately derive H{q±, z) and it is given by 



(9i + 4/) 



(35) 



On the other hand, in the usual collinear factorizationjthe spin-averaged Drell-Yan cross section at small- cc regime is 
dominantly by the qg — > 7*9 channel and is given by [33|, 



da(p^A l+l-X) _ a, 



dx dxA 

X XA 



fq/p{x)fg/A{xA)crqg{s,i,u)S {s + t + U - M^) , 



(36) 



where ctq — 'inal^/SNcsM'^ , s — {P + Pa)^, and (Tqg{s,t,u) has the following form 



d-qg{s,t,u) = 2Tr 



s t 2M''u 
i s si 



with color factor Tr ~ 1/2. Now to connect to the small- a; formalism, realizing 



z{l-z)-' 



t = - 



ql+e^ 



M 



we obtain 



aqg{s,i,u) = 



(37) 



(38) 



(39) 



Then plugging above equation back to Eq. ([55]) . integrating over dxA with the delta function S{s + t + u — M^) 
and using dx/x — dz/z, we immediately end up with an expression exactly the same as the one in Eq. psp. Thus 
we have found that the small-x formalism is exactly the same as the usual collinear factorization formalism for the 
spin-averaged cross section in the relevant kinematic region. 



B. Spin-dependent cross section 



Now let us study whether there is some connection between the spin-dependent cross section dAa/dM'^dyd'^q± in 
the samll-x formalism and that derived in the collinear factorization approach. Again we are interested in the region 
where M^, q^ ^ Ql ^ k\. In the small- a; formalism, it is given by Eq. ([T5)) . For the polarized proton side, 
we have the fc^/j^*(x, fc^), thus we need another linear term from the hard-part functions Ht and Hl to have 
the integral nonvanishing. On the other hand, for the unpolarized nucleus side, to obtain the usual collinear 

gluon distribution function, we need to perform the expansion with respect to just like what we have done for the 
spin-averaged case. In total, we should perform the following expansion, 

dl:^(T{p^A^l+t-X) al^^ ^ 2 ,2, e^'^sj^i 



dAPdyd^qj_ Q-k^AP 



Y^el ^d^k^'^^xf^^\x,kl) J dHd^p^F{xA,p^) 



1 

1 d 



idkldp^^dp 
Then we can perform the d^fc^ integral as follows 



V / d'k^kikyl^^x^kl) = -V / d^k^kl\giytT\^^kl) = i^f r,,^(x,x), (4i) 



where we have used the relation between the quark Sivers function f^'^{x, k'j_) and the twist-3 quark-gluon correlation 
fimction Tq^p{x,x) [l8| 

^ J d^k^klftr^x^kl) = TqAx,x). (42) 
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For the nucleus target side, we perform the (Pp± integral the same way as in Eq. (|32)) for the spin-averaged cross 
section in last subsection. At the end of day, we have 



1 



^ dz 



xTg^F{x,x)xAfg/A{xA)H'^{q±,z), 



where the hard-part function {q±, z) is given by 



8 dk^dpP^dpl 
We could easily work out the explicit expression for H^{q±, z) 

2z^ ( 



[HT{q±,k±,p±,z) + HLiq±,k±,p±,z)] 



fex->-0,Pi-i-0' 



H^iq^,z) = qi 



(43) 



(44) 



(45) 



On the other hand, in the collinear factorization approach, th e sp in-dependent d/S.a / dNP dyd^ q±^ in the region 
-^^j <fi ^ Aqcd is a twist-3 effect and has the following form |33l |34| 



dM^dyd^q^ 



47r2 



^ ' X XA —U 



Tq^F{x,x) - X — Tq^Fix,x) 



mjs, i, u) + Tq,F{x, x)NUs, i it) 



+ Tq^F{x - Xg,x)Hqg{s,i, u) \, 



(46) 



where Xg — —xi/{M'^ — i), both and N^^ are the hard-part functions associated with the so-called soft-gluon pole 
contribution, and -ff^^ is the hard-part function associated with the hard-gluon pole contribution. These hard-part 
functions are given by [H, [sj] 



Hqg{s,i, u) 

N^g{s,i,u) 

H^qisXu) 





s i 2M2-U 




_ i s si 



i£2 



2T, 



R- 



s + u Nf: 



(47) 
(48) 

(49) 



The collinear factorization expression in Eq. (|46p seems very different from the small- a; formalism in Eq. (j43p . 
However, if one considers the so-called forward limit, i.e., \i\ <^ \u\ ~ s, we will find that they are indeed consistent 
with each other as we will show now. From Eq. ([35]), the forward limit is equivalent to q'^ ~ e^j and z — > 1. In this 
limit, the small- a; formalism in Eq. reduces to 



dAajp^A ^ e+e-X) _ „ s 



1 



dM^dyd'^q± 



'^""slq'-^j^-^Y^ el J dzx Tq^F{x, x)xAfg/AixA)- 



(50) 



On the other hand, for the collinear factorization expression in Eq. (|46p . we find that the hard-part function H^g cx 
1/(1 — z) while both N^g and H^g are proportional to 1/(1 — z)^. Thus H^g is sub-leading compared to N^g and H^g 

in the forward limit and can be neglected. At the same time, in the forward limit, Xg — —xi/{M'^ — t) — > thus the 
hard-pole correlation function Tq,F{x — Xg,x) — >■ Tq^F{x,x) is equal to the soft-pole correlation function. Once we 
realize this, we find that the corresponding terms oc N^/{N^ — 1) in Nqg and Hqg cancel each other, only the term 
cx Tr remains. After using the S (^s + i + it — Af^) to take care of the dxA integration and dx/x — dz/ z, we end up 
with the following expression 



dAa{p^A e+l-X) 
dM^dyd^qj_ 



37riV,M2 



2 ^3 9 

q 



dzxTq^F{x,x)xAfg/A{xA), 



(51) 
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which is exactly the same as the forward hmit expression Eq. (j50|) in the smah-a; formahsm. Thus we have demonstrated 
for the spin-dependent cross section that the colhnear twist-three formahsm is consistent with the smah-x formahsm 
in the forward hmit. Since both the spin-averaged and spin-dependent cross sections are consistent with the smaU-x 
formahsm in the forward hmit, the Sivers spin asymmetry Ajsr as a ratio of these two cross sections is thus consistent 
between the small-x formahsm and the coUinear factorization approach. 



V. NUMERICAL RESULTS 



In this section, we present our numerical estimates for the Sivers asymmetry of Drell-Yan lepton pair production in 
both polarized p^-t-p and p''^-|-A collisions at RHIC energies by using Eqs. (fTTl) . (fT8|) . (fT9|) and (|20|) . To be consistent 
with the experimental definition for the asymmetry A^, we choose a frame in which the polarized proton moves 
in the -l-z-direction, the unpolarized proton (or the nucleus) moves in — z-direction, the transverse spin vector s± 
and the transverse momentum of the lepton pair q± along y- and x-directions, respectively. For the transversely 
polarized proton side, we use the quark Sivers function extracted from SIDIS in |36| and flip the sign for the Drell-Yan 
production. For spin-averaged cross section, we use GRV98LO unpolarized parton distribution functions [37| . In 
the current paper, we will not incorporate the evolution effect for the Sivers function [sl, [s^ and will leave such a 
study for a future publication. For the target proton or nucleus side, we need the dipole gluon distribution function 
F{xa,P±), which characterizes the dense gluon distribution in the Drell-Yan process in the proton or nucleus at the 
small-x regime. In the numerical estimate, we choose two parametrizations for the dipole gluon distribution: the 
GBW model in [i^l , and the solution of the BK evolution equation [4l| . 

In the GBW model, the dipole gluon distribution has a Gaussian form [ioj . 



F{XA,P±) = 



1 



where the saturation scale for the proton Q'^{xa) = QIo{xa/xo)_ 
X — 0.288. For the saturation scale of the nucleus, we follow Ref. [42 



^ with QsQ = 
and choose 



;a[xa) =cA^/^Q%{xa/xo) \ 
with c = 0.5 for the minimum bias p+A collisions. 



(52) 

1 GeV, xo = 3.04 x lO"'^, and 



(53) 
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FIG. 3. The single transverse spin asymmetry y4]v is plotted as a function of the transverse momentum q± at rapidity y — 3 
for RHIC energy ^/s = 510 GeV. We have integrated over the invariant mass of the lepton pair 4 < A/ < 8 GeV. The red 
solid curve is for polarized p^-|-A collisions, while the blue dashed curve is for polarized p^+p collisions. On the left plot, we 
have used GBW model for the dipole gluon distribution; while on the right plot, the dipole gluon distribution is taken as the 
numerical solution of the BK evolution equation. 



For the second choice, we use the numerical solution of BK evolution equation for the dipole gluon distribution in 
our calculation. Particularly we take those used in Ref. Q which are able to describe successfully the single inclusive 
hadron production in forward d-|-Au collisions at RHIC in the leading order formalism. Ref. Q has used McLerran- 
Venugopalan (MV) model [43j as the initial condition for the BK evolution equation, which avoids an unphysical 
exponential fall-off of the dipole gluon distribution at large transverse momenta. 
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In Fig. [3l we plot the Sivers asymmetry Ajv of the Dreh-Yan dilepton production as a function of the pair's 
transverse momentum q± at RHIC energy ^/s = 510 GeV. To keep in the forward rapidity region (such that the 
gluon momentum fraction xa in the target is smah thus to justify the use of our smaU-x formahsm), we choose the 
rapidity y — 3. We have integrated the pair's invariant mass M from 4 to 8 GeV. On the left panel, we show the 
result by using the GBW model for the dipole gluon distribution. On the right panel, the asymmetry is calculated 
with the dipole gluon distribution taken as the numerical solution of the BK evolution equation. The blue dashed 
curves are for the asymmetry in polarized p^+p collisions, and the red solid curves are for the asymmetry in polarized 
p^+A collisions. We find that the GBW model and the numerical solution of the BK evolution equations actually 
give very similar results on the Sivers asymmetry for both p^+p and p^^+A collisions. The asymmetry goes to zero 
when the pair's transverse momentum goes to zero q± = 0. This is because that once integrated over both k± and 
p± in Eq. (|18l) , the only remaining transverse vectors arc q± and the spin s± , both of which are needed to form the 
azimuthal spin asymmetry ^ e"^s"g|^. Thus if q± = 0, there are no enough vectors to generate the asymmetry. 
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FIG. 4. The single transverse spin asymmetry ^at is plotted as a function of the rapidity y a,t q± — 1 GeV for RHIC energy 
= 510 GeV. We have integrated over the invariant mass of the lepton pair 4 < M < 8 GeV. The red solid curve is for 
polarized p^+A collisions, while the blue dashed curve is for polarized p^+p collisions. On the left plot, we have used GBW 
model for the dipole gluon distribution; while on the right plot, the dipole gluon distribution is taken as the numerical solution 
of the BK evolution equation. 



We also find that for relatively small q±, the Sivers asymmetry Ajsi becomes smaller from p^+p to p^+A collisions. 
This is expected since the Sivers asymmetry is generated by the small quark transverse momentum k± in the polarized 
proton, whereas the role that the saturation momenta play in this process is to push the produced Drell-Yan lepton 
pair towards larger q± region. Since the saturation scale Qs is larger in the nucleus target, the Sivers effect is 
attenuated and thus the asymmetry is smaller for the nucleus case as compared to the case with proton targets 
for fixed value of q±. On the other hand, when q± becomes relatively large (larger than the saturation scale QsA in 
the nucleus), the saturation effect (due to the difference in Qs) becomes less important, thus we see the asymmetry 
becomes similar for p^+p and p^+A collisions in this region. A comparative measurement of the asymmetry for both 
p^+p and p^+A collisions thus might provide an interesting and novel way to study the saturation scale. 

In Fig.Hl we plot An as a function of the dilepton rapidity y at transverse momentum q± = 1 GeV at RHIC energy 
^/s = 510 GeV. The left plot is for GBW model and the right plot is for the numerical solution of the BK evolution 
equation case. Again y-dependence of the asymmetry is similar for these two cases, and the asymmetry is smaller 
for p^+A collisions than that for p^+p collisions, due to the same reason as in the g_L-dependence case. It is easy to 
understand the reason why the asymmetry is peaked at y ~ 3, which corresponds to the parton momentum fraction 
in the polarized proton Xp 0.2. This effect is simply due to the fact that the parametrization we used for the Sivers 
function has a maximum at x ~ 0.2 [36|. Nevertheless, we also find that the saturation effect plays a non-trivial role 
in the calculation of the asymmetry. 



VI. SUMMARY 



We studied the Sivers single spin asymmetry of the Drell-Yan lepton pair production in the forward rapidity region 
(small-x regime) for both polarized p^-l-p and p^-|-A collisions. By including the transverse momentum dependence of 
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the quark distribution inside the polarized proton, we calculated both the spin-averaged and spin-dependent Drell-Yan 
differential cross sections in the small- a; formahsm. In the so-called geometric scaling region in which the invariant mass 
M of the pair is much larger than the saturation scale Qs and the pair's transverse momentum q±, we demonstrate 
that the Drell-Yan cross section in the small-x formalism agrees completely with those derived from the transverse 
momentum dependent factorization approach. We further show that the small- a; formalism is also consistent with 
the coUinear factorization formalism in the so-called forward limit. We evaluated the SSAs of Drell-Yan lepton pair 
production in both polarized p^-l-p and p^-|-A collisions for RHIC energy ^/s = 510 GeV at forward rapidity region. 
We find that in the region when the pair's transverse momentum q± is small, the Sivers asymmetry is smaller in p^-|-A 
collisions as compared to p^-l-p collisions, while the asymmetry becomes similar at relatively large q±. The future 
measurement and comparative study of the Sivers single spin asymmetry for the Drell-Yan lepton pair production in 
both polarized p^-l-p and p^-|-A collisions will be very interesting and important as they probe both the transverse 
spin and the small- a; saturation physics. 
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